Resonant dynamics of H atom in elliptically polarized microwave field 
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The dynamics of Rydberg states of atomic hydrogen driven 
by elliptically polarized microwaves of frequency fulfilling 
2:1 classical resonance condition is investigated both semi- 
classically and quantum mechanically in a simplified two- 
dimensional model of an atom. Semiclassical results for 
quasienergies of the system are shown to be in a good agree- 
ment with exact quantum data. The structures of the quan- 
tum states are found to reflect the underlying classical dy- 
namics; especially we show the existence of nonspreading 
wavepackets propagating on elliptical trajectories. 

PACS: 05.45. +b, 32.80.Rm, 42.50.Hz 

The pioneering experiment on Rydberg hydrogen (H) 
atoms ionization in microwave fields [[l] and a subsequent 
interpretation of its results in terms of the chaotic clas- 
sical dynamics opened up a new possibility to study 
dynamics of atomic systems whose classical counterpart 
reflects a chaotic behavior (for review see ||). Quantum 
mechanically the behavior of atoms in microwave fields 
is associated with multiphoton processes what makes 
the quantum perturbation calculations hardly possible 
(many atomic states are strongly coupled) while, on the 
other hand, allows to apply semiclassical or even purely 
classical description. Signatures of local structures in 
classical phase space have been observed in the experi- 
mental ionization thresholds where an existence of stable 
islands caused increase of the stability of atoms for ion- 
ization S. 

The stable structures are suitable for a semiclassical 
quantization. The largest of them are created by the pri- 
mary resonances between the driving field and the un- 
perturbed Kepler motions of the Rydberg electron. In 
terms of the principal quantum number no of the un- 
perturbed atomic state, the ratio of the microwave fre- 
quency, u>, to the Kepler frequency, lok = V n o> nas to 
fulfill lu/lok = s, where s is an integer number. Some 
of states supported by the resonance islands reveal inter- 
esting properties from the quantum-classical correspon- 
dence principle point of view. Their time-dynamics may 
be viewed as a motion of non-dispersive wavepackets as 
shown for both circular polarization (CP) [0 and linear 
polarization (LP) || of microwaves. These new states 
attracted a lot of attention, their decay via ionization || 
or spontaneous emission as well as ways to populate 
such states in experiments || have been studied. The 
ways of modifying wavepackets properties by additional 



static fields were also considered p,|l0[. Semiclassically, 
all of them are supported by the principal resonance is- 
land, i.e. uj/luk — 1. The wavepackets built inside a 
higher primary resonance island have been discussed by 
Holthaus in model one dimensional (ID) systems and 
to a realistic H atom in LP microwaves |l2| . 

Theoretical detection of such interesting states obvi- 
ously needs reliable semiclassical methods. In the LP 
case the ID approximation to the dynamics allowed to 
quantize resonant states exploring Mathieu equation jl3) . 
In the CP case passing to the rotating frame removes ex- 
plicit time dependence of the system and the harmonic 
approximation around the stable fixed point gives the ap- 
propriate semiclassical predictions . For the realistic 
LP problem as well as for arbitrary elliptical polarization 
(EP) in two dimensional (2D) model of an atom the semi- 
classical method based on the Born-Oppcnhcimer ap- 
proximation allows to describe resonant dynamics • 

This brief report extends the previous analysis of the 
EP case [[u], limited to the principal resonance island, 
to a higher primary resonance. We describe the full dy- 
namics of quasienergies as a function of the microwave el- 
lipticity for the 2:1 resonance case, concentrating in par- 
ticular on nonspreading wavepackets. The experimental 
data for the EP problem are available in the range of the 
scaled frequency, luq = w/wx, up to 1.4 Thus a the- 
oretical description of the system for ujq — 2 should be 
able to be compared with an experiment probably in the 
immediate future. 

As previously jL4|, we treat semiclassically and quan- 
tum mechanically a simplified 2D model of an atom. 
Study of such simplified models have been most success- 
ful in the past both for the LP problem (where ID model 
has been a main source of quantum results for a long 
time ||) and for the CP one where also the 2D, polariza- 
tion plane restricted model has been utilized (and 
references therein). 

The Hamiltonian of the hydrogen 2D model atom 
driven by an elliptically polarized electromagnetic field 
reads in the dipole approximation (in atomic units) 

p 2 + p 2 1 

H=— - h Ftxcosuit + ays'mujt), (1) 

2 r 

where r = \/ x 2 + y 2 while F and cu denote the amplitude 
and the frequency of the microwave field, respectively. 
a defines the ellipticity of the microwaves, with a = 
(a = 1) corresponding to a LP (CP) limiting case. 
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Using the Floquet theorem pq] the solution of the 
quantum problem is equivalent to diagonalizing the Flo- 
quet Hamiltonian, TL = H — id/dt with time-periodic 
boundary conditions, getting the eigenvalues (quasiener- 
gies) and time periodic eigenstates (Floquet states) . The 
semiclassical quantization of resonant dynamics, based 
on a prescription of [Til , closely resembles the similar 
procedure applied by us recently for the LP and EP case 
p2|p^ . The method bases on passing to the extended 
phase space by defining the momentum p t conjugate to 
the t (time) variable what yields the new Hamiltonian, 
TL = H + p t (l^] . The quasi-energies of the system will 
be then the quantized values of TL. As the next step we 
express the Hamiltonian in action-angle variables of the 
unperturbed Coulomb problem jl8). For the 2D model 
atom those are e.g. the canonically conjugate pairs (J, 9) 
and (L, cf>). J is the principal action (corresponding to the 
principal quantum number, no). The conjugate angle, 9, 
determines the position of the electron on its elliptic tra- 
jectory and depends linearly on time, 9 ~ toxt, for an 
unperturbed atom. L is the angular momentum (equal 
to L z for the 2D motion in the x — y plane) while <fi is the 
conjugate angle (the angle between the Runge-Lenz vec- 
tor and the x-axis, i.e. the main axis of the polarization 
ellipse). 

Considering the case of the resonant driving, i.e. luq — 
s, we apply the secular perturbation theory [[l7j to aver- 
age over the nonresonant terms which yields the approx- 
imate resonant Hamiltonian of the form 



TL r 



where 
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2s 2 J 2 



loJ + FT(L, (f>; a) cos[6» - f3{L, 0; a)} + p t 

(2) 



= s9 — tot, J = — , p t 
s 



-■ Pt + <jjj- 



(3) 



The explicit form of T(L, (f>; a) and f3(L, a) is given by 
Eq. (2.17) of @. 

The last stage is to quantize the system using the ap- 
proximate Hamiltonian, Eq. (^). Trivial quantization of 
p t , exploring the time periodicity the system (note that 
in the rotating frame, defined by Eq. (|^), the time period 
is r = 2sir/oj |fl2"||), yields additive terms kui/s to quasi- 
energies, where k is an integer number [lljjl^] . Thus the 
spectrum associated with states localized in the s : 1 
resonance island repeats itself along the energy axis at 
distances lu/s [[ll 12\. The radial motion in ( J, 9) space 
is much faster than the angular motion in (L, (j)) space 
[^9|,Q. Hence, in the spirit of the Born-Oppenheimer 
approximation, one may first quantize the fast radial mo- 
tion, keeping L and (j) fixed, then pass to the quantization 
of the slow angular motion. In fact, because of the spe- 
cific form of Eq. , the order of the quantizations does 
not matter, one may quantize first the slow motion yield- 
ing discrete values of T and then go to quantization of 



the fast motion |L^ , |T4[| , We consider below the s = 2 
resonance as a generic example. 

We are interested in strongly localized, wavepacket-like 
states lying close to the center of the resonance island 
in fast (J, 9) variables. Previously we have used 

pure WKB quantization for that motion which, however, 
works poorly when the island size is small. Such is a 
case for s > 1 resonances. Thus we improve the pro- 
cedure and expand the principal action to the second 
order around the center of the resonance island defining 
I = J — (s 2 ^) -1 / 3 . Eq.(||) gives then a standard pendu- 
lum Hamiltonian with the island size given by \J FT. To 
ensure a maximal radial localization of the electron we 
consider then a ground state of the fast motion by tak- 
ing appropriate quantum eigenvalues (as given by Math- 
ieu equation solutions |20[ ), see similar treatment for ID 
systems 

The slow angular motion is determined by constant 
values of T(L, (f>; a). To compare the semiclassical pre- 
dictions to the exact quantum calculations, we consider 
the no = 42 manifold of our 2D model atom. For res- 
onant driving we take the microwave frequency to be 
to = 2uj k = 2/(n + 1/2) 3 , i.e. uj = 2. Note that in 
the 2D model the effective principal quantum number is 
half integer. Fig. [j] shows values of T as a function of 
the scaled angular momentum, La — L/(no + 1/2), and 
the (j> angle for two different values of the field elliptic- 
ity, a. Semiclassically quantized contours, reflecting slow 
evolution of the electronic ellipse, are also presented in 
the figure. Those of them which are localized around 
extrema of T correspond to states with well defined elec- 
tronic ellipse. Size of the resonance island in the (J, 9) 
space depends on value of T, thus, only states localized 
around the maxima will show strong radial localization 
too. Note that T is equal to zero for circular orbits, i.e. 
for Lq = ±1. It is obvious because circular motion is 
purely harmonic and no primary resonance exists except 
the 1 : 1 case. 

Consider level dynamics with a change of the field el- 
lipticity. Fig. || shows semiclassical and numerical results 
for quasienergies corresponding to the resonantly driven 
n = 42 manifold as a function of a, for the scaled field 
amplitude Fq = F(no + 1/2) 4 = 0.03. The nice quan- 
titative agreement between semiclassics and numerics is 
achieved (except in the region of broad avoided crossings 
with other levels - partners in the crossing are not plot- 
ted for clarity) with significant improvement over earlier 
approach |L4| for small T region. 

This level dynamics is easy to understand by inspec- 
tion of angular motion change with a, Fig. ^. Semiclas- 
sically, for a = 0, all states are degenerate because of the 
symmetry of the (L, <fi) space with respect to the L = 
axis. The highest degeneracy exists for librational states 
situated in elliptical islands around Lq ~ ±0.6, = 0,7r 
because the islands are identical. With a small change 
of a, values of T corresponding to negative Lq become 
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smaller while those corresponding to positive Lq become 
greater. Quasienergy levels simply follow the increase or 
decrease of T, i.e. the greater value of T the higher the 
corresponding quasienergy level. Thus, for a > 0, the 
degeneracy of many states is removed. Still for a < 0.17 
there exist three pairs of identical elliptical islands which 
support identical semiclassical states, see Fig. [l]. 

With further increase of a the islands situated around 
Lq rs —0.6, cf> = 0, 7r shrink and finally disappear — they 
support fewer and fewer librational states which during 
an increase of the field ellipticity vault over separatrix 
and become rotational ones. Also the islands situated, 
for a = 0, around Lq — 0, <j) — 7r/2,37r/2 shrink with a 
change of the ellipticity. Additionally they move towards 
higher negative values of Lq. For a close to unity all 
elliptical islands are too small to support semiclassical 
states. Then all states are rotational. 

As mentioned before strong localization in both angu- 
lar and radial motion, supporting the existence of non- 
spreading wavepackets is expected around maxima of T. 
Thus one expect nonspreading wavepacket character for 
states localized around Lq w ±0.6, <p = 0, ir for a < 0.17. 
For the greater field ellipticity only the islands around 
Lq w 0.6, 4> — 0, 7r could support nonspreading wavepack- 
ets. However, single wavepacket propagating along the 
2:1 resonance periodic orbit could not fulfill the period- 
icity of Floquet states because the period of the orbit 
is twice longer than the microwave period. So one ex- 
pects Floquet states being linear combinations of two 
wavepackets shifted in t, i.e. a coordinate variable in 
the extended phase space, by 2ir/tu which exchange their 
positions after the microwave period |ll],[l2| • 

As a representative of such wavepackets propagating 
on an elliptical trajectory we have chosen to plot the 
states localized around Lq w —0.6, <f> — 0, ir for a = 0.1. 
Because of the tunneling effect a single quantum eigen- 
state contains a symmetric or antisymmetric combina- 
tion of two semiclassical solutions corresponding to the 
ellipses with the Runge-Lenz vector directed parallel or 
antiparallel to the x-axis. Linear combinations of two 
such eigenstates allow to remove one of the ellipses. The 
resulting state is localized on a single ellipse but still it 
consists of symmetric combination of two wavepackets 
shifted in microwave phase by 27r (for the 2:1 resonance 
case To separate a single wavepacket one has to 

find another state localized on the same ellipse but con- 
taining antisymmetric combination of the wavepackets. 
The desired state is prepared using eigenstates coming 
from the similar manifold shifted by w/2 jll],[l2]]. The 
resulting wavepacket moving on an elliptical trajectory 
is shown in Fig. |^. This is the single wavepacket rotat- 
ing in the opposite direction to the direction of the field 
vector rotation. It propagates along the periodic orbit 
supported by the 2:1 resonance, thus, the period of the 
motion is twice longer than the microwave period. 

In conclusion we would like to stress that the 



wavepacket presented in Fig. |3| is not an eigenstate of the 
system. Tunneling effects between the stable classical is- 
lands will be changing the shape of the packet but, at 
least in the example considered, the corresponding time 
scale is of order of a few hundreds of the microwave pe- 
riod. Another mechanism of its destruction is a slow 
ionization 

The analysis presented is restricted to the 2D model, 
its validity for the real three dimensional atom is an open 
question. Certainly, in the limiting LP case, due to the 
azimuthal symmetry, the wavepackets appear as dough- 
nuts shaped localized functions moving up and down (as- 
suming a vertical polarization of LP microwaves) |p.2| . 
For the CP case, on the other hand, the wavepacket mo- 
tion was found to remain essentially 2D fl. The inter- 
esting problem how the third dimension affects the dy- 
namics for the general EP case is left for future work. 

Finally we would like to note that the 2:1 resonance is 
an example of a general s : 1 resonance as the angular 
motion generated by T is topologically the same for s > 2. 
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FIG. 1. Two dimensional hydrogen atom illuminated by 
resonant, u>o = 2, elliptically polarized microwaves. The ef- 
fective scaled perturbation To = Ylo 2 ^ 3 is plotted as a function 
of the scaled angular momentum, Lq, and the angle, (j>, be- 
tween the Runge-Lenz vector and the main axis of the polar- 
ization ellipse - upper row. Bottom row shows equipotential 
curves of the angular part V of the Hamiltonian Ti r , Eq. (jij), 
representing the slow evolution of the Kepler ellipse. The 
curves correspond to semiclassical states originating from the 
no = 42 hydrogenic manifold. Columns correspond to the dif- 
ferent ellipticity of the microwaves a = 0.1 (left) and a — 0.5 
(right). 

FIG. 2. Two dimensional hydrogen atom driven by reso- 
nant, wo = 2, elliptically polarized microwaves. Level dy- 
namics, versus a (i.e. the degree of the field ellipticity), of 
the semiclassical quasi-energies [panel (a)] of the states orig- 
inating from the no = 42 hydrogenic manifold for Fo = 0.03 
compared with the exact quantum results [panel (b)]. 

FIG. 3. Wavepacket, being a linear combination of 4 eigen- 
states of the hydrogen atom plus elliptically polarized mi- 
crowaves system with the field amplitude Fo = 0.03, fre- 
quency ujo = 2 for no = 42 and the ellipticity a = 0.1. Tem- 
poral evolution is plotted at times cut = (top left), it/2 (top 
center), it (top right), 37r/2 (bottom left), 2tt (bottom right). 
This wavepacket rotates on an elliptical orbit in the opposite 
direction to the rotation of the field vector and essentially 
repeats its periodic motion with period 4tt/uj. It slowly dis- 
perses, either because the 4 building states are not exactly 
degenerate (tunneling effect) or because it ionizes. The size 
of each box is ±4000 Bohr radii in both x and y directions. 
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